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Abstract. We examine the momentum dependence of the bispectrum of two-field inflationary 
models within the long-wavelength formalism. We determine the sources of scale dependence 
in the expression for the parameter of non-Gaussianity /nl and study two types of variation 
of the momentum triangle: changing its size and changing its shape. We introduce two 
spectral indices that quantify the possible types of momentum dependence of the local type 
/nl and illustrate our results with examples. 
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1 Introduction 

Non-Gaussian features of the cosmic microwave background have been an important subject 
of cosmological research in recent years. On the one hand non-Gaussianity has not yet been 
observed after seven years of WMAP [1] and thus is expected to be small, on the other hand 
the PLANCK satellite [2] is expected to give precision constraints, and possibly detect it, 
soon. The importance of a discovery or further constraints lies in the fact that it will enable 
us to discriminate between different classes of models of inflation since these produce different 
types and amounts of non-Gaussianity. 

There are basically two distinct types of non-Gaussianity that are most important from 
the point of view of inflation: the equilateral type produced at horizon-crossing, which has 
a quantum origin and is maximal for equilateral triangle configurations [3], and the local 
type produced outside the inflationary horizon due to the existence of interacting fields. The 
latter is maximized for squeezed triangles, i.e. isosceles triangles with one side much smaller 
than the other two [4, 5]. The first type is known to be slow-roll suppressed for single-field 
models with standard kinetic terms and trivial field metric [6]. On the other hand, models 
with non-standard kinetic terms coming from higher-dimensional cosmological models are 
known to produce non-Gaussianity of the equilateral type so large that it is not compatible 
with WMAP observations [7-10], thus leading people to consider an extra field in order to 
achieve smaller values of the parameter /nl of non-Gaussianity. 

Non-Gaussianity of the squeezed type can be found naturally in multiple-field models of 
inflation [11, 12], due to the sourcing of the adiabatic mode by the isocurvature components 
outside the horizon. For single-field models this is obviously impossible due to the absence 
of isocurvature modes. There has been much study of two-field models [13-20], being the 
easiest to investigate, in the hope of finding a field potential that can produce local non- 
Gaussianity large enough to be measurable. It proves to be non-trivial to sustain the large 
non-Gaussianity produced during the turn of the fields until the end of inflation. 
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Non-Gaussianity produced at horizon-crossing is known to be momentum-dependent. 
The scale dependence of the equilateral /nl produced for example from DBI inflation [9, 21- 
24], has been examined both theoretically [25-28] and in terms of observational forecasts 
[29, 30]. In this paper we are going to study the scale dependence of local-type models 
that has not been studied as much. Squeezed-type non-Gaussianity, produced outside the 
horizon, is usually associated with a parameter of non-Gaussianity /nl that is local in real 
space, and therefore free of any explicit momentum dependence, defined through C( x ) = 
Cl(x) + (3/5)/^ ) L aI (CL( a; ) 2 — (Cl{x)) 2 ), where Cl is the linear Gaussian part. Nevertheless, 
calculations of /nl for several types of multiple-field models (see e.g. [19, 20, 31]) show 
that there is always a momentum dependence inherited from the horizon-crossing era, which 
can in principle result in a tilt of /nl- When a physical quantity exhibits such a tilt one 
usually introduces a spectral index, as for example in the case of the power spectrum. The 
observational prospects of the detection of this type of scale dependence of local /nl were 
studied in [30]. Only recently spectral indices for /nl were defined in [32-34], keeping constant 
the shape of the triangle or two of its sides, within the 5N formalism. Note, however, that 
most theoretical predictions have considered equilateral triangles for simplicity, even though 
the local-type configuration is maximal on squeezed triangles. If one were to calculate a 
really squeezed triangle, then /^ al acquires some intrinsic momentum dependence due to 
the different relevant scales, as was shown in [20]. 

It is both these effects we want to study in this paper: on the one hand the tilt of /nl 
due to the background evolution at horizon-crossing and on the other hand the impact of 
the shape of the triangle on /nl- In order to do that in a concrete way, such that these 
effects do not mix, we define two independent spectral indices, each one quantifying different 
deformations of the momentum triangle. Moreover, having an exact expression of /nl for 
an isosceles triangle, we are able to study and understand for the first time the origin of 
both types of momentum dependence of /nl- We also provide analytical estimates for the 
quadratic model (which actually hold for any equal-power sum model) that we use in this 
paper to illustrate our results. 

The paper is organised as follows. We begin by reviewing the power spectrum and the 
bispectrum along with the background of the theory in section 2. In section 3 we present 
the long-wavelength formalism results and discuss the sources of scale dependence in the 
expression for /nl- We also introduce two spectral indices, able to quantify the effects of 
different triangle deformations. In section 4 we study the scale dependence for triangles of 
constant shape but of varying size, which is mainly due to horizon-crossing quantities, while 
in section 5 we study the scale dependence related to the shape of the triangle. Finally, we 
conclude in section 6, while analytical expressions for the final value of the spectral indices 
for any equal-power sum potential are given in Appendix A. 

2 Preliminaries 

The power spectrum of the cosmological adiabatic perturbation £i has proved to be a valuable 
tool to connect the theory of inflation to observations of the sky today. The power spectrum 
provides two important observables: its amplitude V^{k,t) 



carrying information about the mass scale of the inflaton, and the spectral index quanti- 
fying the tilt of this almost constant, scale-independent spectrum due to the small (assuming 




(2.1) 
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slow roll) but non-zero evolution of the inflationary background, 

n _ i = dlnV < = dlnV < dt * = dlnV ( 1 (2.2) 
d In k dtk d In k dtk 1 — e& 

From now on a A; subscript on the slow-roll parameters will denote evaluation at the time 
that the scale k exits the horizon. The last equality in (2.2) comes from the relation k = aH 
at horizon crossing, while using 

H = -eH, 

a = a. (2.3) 

Here we choose the time coordinate to be the number of e- foldings (for details see [20]). The 
first of the above equations is just the definition of the e slow-roll parameter, while the second 
is the form that the definition of the Hubble parameter acquires due to the choice of the time 
variable. 

We also give here the Einstein and field equations [35] for the inflaton fields 4> A , with 
the index A numbering the fields, that roll down a potential W: 

(2.4) 





k 2 /n 2 


H 2 = 






3 V 2 


tl A = 


-m A - 



H 



with K 2 = 8irG = Siv/m^ and W t A = dW/d<j> . T1 A = H<j) A is the canonical momentum of 
the inflaton. Here we choose to work with two scalar fields, so that A = 1,2, with canonical 
kinetic terms and a trivial field metric. Generalizing any of these assumptions would still 
allow for the system to be studied in principle, although the computations would be more 
complicated. 

Apart from e, a hierarchy of slow-roll parameters can be built from the time derivatives 
of the inflation fields: 

rfn)A _ _J_ ( h^Y' 1 U A (2 . 5) 

We will denote the first and second-order slow-roll parameters of the hierarchy as r/^ A = rj A 
and rj^ A = £ A . In order to distinguish adiabatic from isocurvature effects we construct a 
basis in field space {e A ,e A }, where e A is a unit vector parallel to the velocity of the fields 
and e A a unit vector parallel to the component of the acceleration that is perpendicular to 
the velocity: 

ef = —— and e, = -=— ; -. : — . (2.6) 

1 n 2 \U A -e A e 1B U B \ V ^ 

This is why we use the 1 subscript for the adiabatic perturbation £i, while the isocurvature 
perturbation will be denoted as £2- The projection of the slow-roll parameters on e A will be 
denoted by || and the projection on e A by _L. Finally we also define the slow-roll parameter 

X as 

W22 . _!_ II (O 7\ 

where W22 = W' AB e2A^2B is the second derivative of the potential W of the fields projected 
in the 22 direction. 
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In addition to the power spectrum we can gain more information from the CMB by 
studying the Fourier transform of the three-point correlation function, 

(Ci fcl Cifc 2 Cifc 3 > = {2tt)^I 2 5\Y, k s )B c (h,k 2 , k 3 ), (2.8) 

s 

where is the bispectrum. Because of the overall J-function we see that the vectorial sum 
of the three k- vectors has to be zero. In other words, the three A;- vectors form a triangle. The 
amplitude of the bispectrum can provide additional constraints on the slow-roll parameters 
of a given type of inflationary model. The profile of the bispectrum, i.e. the shape of the 
momentum triangle, gives information on the type of the inflationary model itself. For 
example, models with higher-order kinetic terms produce a bispectrum of the equilateral 
type (see e.g. [36]), mainly due to quantum interactions at horizon crossing. By equilateral 
type we mean here a bispectrum that becomes maximal for equilateral triangles. On the 
other hand, canonical multiple-field inflation models predict a bispectrum of the local type. 
This arises from non-linearities of the form £1 = £il — (3/5)/nl(Cil ~~ (Cil) 2 ) (Cil being the 
first-order adiabatic perturbation) that are created classically outside the horizon, leading to 
a bispectrum of the form 



B ( (k u k 2 ,k 3 ) = -jj/ NL (^-'P c (k 1 )^-P c (k 2 ) + (k 2 o fc 3 ) + (A* o fafj , 



(2.9) 



where /nl is usually assumed to be constant. This bispectrum becomes maximal for a 
squeezed triangle, i.e. a triangle with two sides almost equal and much larger than the third 
one. As we will discuss in the rest of the paper, /nl is not actually a constant, but depends 
on the size and shape of the momentum triangle. 

In order to study the dependence of the non-Gaussianity on the shape of the triangle, 
instead of using ki,k 2 , and k 3 we will use the variables introduced in [37, 38], 

T s h + k 2 + k 3 h-k 2 k 3 - ki - k 2 

K = - , 7=-^, P = 2# ' (2 - 10) 

which correspond to the perimeter of the triangle and two scale ratios describing effectively 
the angles of the triangle. They have the following domains, < K < oo, < (3 < 1 and 
— (1 — /3) < 7 < 1— (3, see figure 1. As one can check from the above equations, the bispectrum 
becomes maximal for (3 = 1 and 7 = 0, or (3 = and 7 = ±1, i.e. for a squeezed triangle. 
In this paper we always assume k\ = k 2 , dealing only with equilateral or isosceles triangles 
(note that the relation k\ = k 2 is satisfied by definition for both equilateral and squeezed 
triangles). The two scales of the triangle k 3 = k<k' = k\ = k 2 can be expressed in terms of 
the new parameters f3 and K as 

k = {l-p)K and k' = ^-^-K. (2.11) 

In this paper we are going to study only acute isosceles triangles 1/3 < (3 < 1, but the 
generalization is straightforward. 

3 Sources of scale dependence 
3.1 Long- wavelength results 

In this paper we use the long- wavelength formalism to study the parameter of non-Gaussianity 
/nl and its scale dependence. The long-wavelength formalism consists of doing a perturbative 
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Figure 1. The types of deformations of the momentum triangle we are considering. Left: Conformal 
transformation of the triangle. In the 7, j3 plane this corresponds to a constant (red) point. Right: 
Keeping the perimeter of the triangle K — const we change the shape of an isosceles 7 = triangle, 
moving along the bold red line for k± = > &3- 



expansion of the exact differential equations for the non-linear cosmological perturbations on 
super-horizon scales, which were obtained by neglecting second-order spatial gradient terms 
(see [11, 20, 35]). Slow-roll solutions for the first-order perturbations at horizon-crossing are 
used as initial conditions. 1 The super-horizon assumption is equivalent to the leading order 
of the spatial gradient expansion and requires the slow-roll assumption to be satisfied at 
horizon-crossing (but not afterwards, see [20, 40, 41]). The non-Gaussianity parameter for 
an isosceles triangle of the form k\ = k 2 = k' > £3 = k was found in [20] to be 



6, -2v 12k ,/[l + {vi2k'f 



NL 



Ik 



vi2k'\9sr(k ,k )+g iso (k ,k )+g in t(k ,k) 



+2%fi2fc 9sr(k', k)+g iso (k', k)+g int (k', k) 



, (3.1) 



where /nl = /nl(*; tk', tk) depends on t k i and t k , denoting the horizon-crossing times of the 
two scales k and k! of the triangle, respectively. This result is exact and valid beyond the 
slow-roll approximation after horizon-crossing. All quantities appearing in this formula will 
be explained below. 

The quantity v\ 2 is a transfer function showing how the isocurvature mode (denoted 
by the subscript 2) sources the adiabatic component C,\. In the following two more transfer 
functions appear, namely v 22 and €>32 showing how the isocurvature mode sources the isocur- 
vature component Q 2 and the velocity of the isocurvature component 9 2 = C, 2 , respectively. 
v a2 is a function of the horizon-exit time t k of the relevant perturbation of scale k and it 
also evolves with time t, at least during inflation. In (3.1) as well as in the formulas that 
follow, v a2 k = v a2 (t,tk)- The indices a, b take the values 1,2,3, indicating respectively the 
adiabatic perturbation (j, the isocurvature perturbation £21 and the isocurvature velocity 
&2- 2 v a2 comes from the combination of the Green's functions G a2 and G a z of the system of 



1 The non-Gaussianity produced at horizon-crossing is not included in the long- wavelength formalism, but 
can be computed in another way, see [39], and added by hand. It is negligibly small in models with standard 
kinetic terms. 

2 Due to the exact relation Q\ = 2t7 X ^2, there is no need to consider the velocity of the adiabatic perturbation 
9i as an additional variable [11]. 
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equations for the super-horizon perturbations (for the system of equations that the Green's 
functions obey see [11, 20]): 

MMk) = G a2 (t,t k ) - X kG a3 {t,t k ). (3.2) 

The quantity ^ k in (3.1) is defined as 7^ = —KH k /(2k 3 ^ 2 y /e k ). 

Except for the overall factor, /nl has been split into three contributions: g sr , gi so 
and gint- 3 g sr is a term that is slow-roll suppressed, since it depends only on horizon-exit 
quantities, where by assumption slow-roll holds, 

n 7 \ ±( G 22klk2 vi 2kl 1 G 22klk2 \ , 3 Xfc 2 n 3 || 

V 2 7J 12 fc 2 ^^12fei J 4 2 1 

H — 7- *z h <-J22fcifc 2 — 

4 ^ v 12k2 J 2(v 12 ) 2 

. G\ 3 (t,t kl ) 



2 



3 {Xk±G 22klk2 - Xk 2 G 33kl k 2 ) . ri f 3 + e fci+ 2? ?fc 1 . J_ 

T= + CT32fc lfc2 — — + rj kl 

2fl2fci \ ^12fci 



~^G 32klk2 - lrG 12klk2 I e kl + 77J + 2^ - (1 + + *- ^ J • (3-3) 

4 ^ y * vi 2kl j 

Here we introduce some new notation, 

(v 12 ) 2 = vi 2kl v uk2 , (v 22 ) 2 = v 22kl v 22k2 , (v 32 ) 2 = v 32kl v 32k2 , 
v 22 v 32 = \{v 22kl v 32k2 + v 22k2 v 32kl ), 

and also G abklk2 = G a b(t kl ,t k2 ). Moreover, we assume k\ > k 2 . g sr is the only term from 
which a (small) part survives in the single-field limit, i.e. in the limit where vi 2 = at all 
times. For the equilateral case k' = k the two last lines of g sr are zero, since the Green's 
functions satisfy 4 

G ab (t,t) = 5 ab . (3.4) 
The contribution gi so is a term that survives as long as the isocurvature modes are alive, 

9iso{h,k 2 ) = (e + ??")(^22) 2 + v 22 v 32 . (3.5) 

If at the end of inflation these are non-zero, /nl can still evolve afterwards and we cannot 
be sure that its value survives until today. Finally, gi n t is given by 

9int(h,k 2 )=- [ dt'\2{i 1 ± ) 2 {v 22 ) 2 +{e+r,\\)v 22 v 32 +^ (3.6) 
Jt kl 1 J 

with 

C = U^x - 6r/V + 6(7/11) V + G^) 3 - 277^11 - 2r ? "e ± - J\-^p (W211 + W 222 ), (3.7) 



3 In [20] we had also a fourth contribution gk, denoting the terms that vanish for an equilateral triangle. 
Here we have incorporated these terms in g 3r (the last two lines), since they are also slow-roll suppressed. 

4 In addition, in the equilateral case k' = k the 7^ ratios in (3.1) reduce to 1 and the two terms in the 
brackets of (3.1) become identical (apart from the factor 2). 
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Figure 2. The relative change of the horizon-crossing first-order slow- roll parameters e (bold curve), 
77" (dot-dashed curve), r/ 1 - (dotted curve) and x (dashed curve) at tk as a function of the ratio k/kso 
of the horizon-exit scale to the scale that left the horizon 60 e-folds before the end of inflation, for the 
model (3.9) with mass ratio m^/m a = 9. 

where W mn i = W' ABC e m Ae n Beic- It is from this integrated effect that any large, persistent 
non-Gaussianity originates, if we consider only models where the isocurvature modes have 
vanished by the end of inflation. For the analytical approximations that we will provide 
(in addition to the exact numerical results), it is useful to note that within the slow-roll 
approximation gi n t can be rewritten as 



where gi n t is another integral that is identically zero for the two-field quadratic model, or 
even more generally for any two-field equal-power sum model (see [20] for details). 

3.2 Discussion 

Inspecting (3.1) ones sees that there are two sources of momentum dependence for /nl : the 
slow-roll parameters at horizon-crossing and the Green's functions G a i, or their combinations 
v a 2- In order to study their impact we shall use the quadratic model 



with m^/rriu = 9. The procedure to follow is to solve for the background quantities (2.4) and 
then for the Green's functions (see subsection 2.2 in [20]) in order to apply the formalism. The 
quadratic model's Green's functions can be found numerically, or even analytically within 
the slow-roll approximation, which is valid for a small mass ratio like the one we chose here. 
However, all our calculations in this paper are numerical and exact, without assuming the 
slow-roll approximation after horizon crossing. We only use the slow-roll approximation 
after horizon crossing for the analytical approximations that we provide (e.g. eq. (4.1)) and 
sometimes to clarify the physical interpretation of results (e.g. the use of (3.11) below to 
explain the behaviour of ^12). Inflation ends at tf defined as the time when e/ = 1. From 



9int{kl,k 2 ) = Vi2k 1 G 2 2k 1 k 2 -Vki + 





(3.9) 
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Figure 3. The evolution of the transfer functions v\2 (left) and V22 (right) as a function of the number 
of e-foldings N — Nqo for the time interval around the turning of the fields and for different horizon 
exit scales, varying from top to bottom as kgo x 10, £^0, ^60/10, for the model (3.9) with mass ratio 
m^/m a = 9. 



now on a subscript / will denote quantities evaluated at the end of inflation. We also define 
the scale that exited the horizon 60 e-folds before the end of inflation as k^ and use it 
as a reference scale, around which we perform our computations (fcgo being the scale that 
corresponds to the text books' minimal necessary amount of inflation). 

In figure 2 we plot the first-order slow-roll parameters for a range of horizon-crossing 
times around keo- While the heavy field rolls down its potential the slow-roll parameters 
increase, reflecting the evolution of the background. This implies that /nl,j which is in 
general proportional to the slow-roll parameters evaluated at tk and ty, should increase as 
a function of k and k! . This can easily be verified for the initial value of /NL,m at t = ty, 
which according to (3.1) with v±2k' = takes the value 

£fc' + %, + Vk'G22k'k- (3.10) 

1 + 2 i[i + (GW fc ) 2 ] 

Apart from the slow-roll parameters the other source of momentum dependence for /nl 
lies in the Green's functions and particularly how their time evolution depends on the relevant 
horizon-crossing scale. The two main quantities that we need to study in order to understand 
their impact on /nl are the transfer functions v\2 and r>22- This is due to the fact that U32 is 
slow-roll suppressed and the rest of the Green's functions appearing in (3.1) can be rewritten 
in terms of v\2 and V22 within the slow-roll approximation (for details, see [20]). In particular 
G a 3 = G a2 /3, G 32 (t,ifc) = -x(t)G22(t,tk) and hence G a2 « v a2 . Note that except for the era 
of the turning of the fields, the slow-roll assumption is a good approximation during inflation 
in this particular model. The slow-roll evolution equations for vi2k and V22k are 

^vi2k = 2r? ± 'U 2 2fc and ~j^ 22fc = ~X^22k- (3.11) 

As was discussed above, v±2 describes how the isocurvature mode sources the adiabatic 
one, while V22 describes how the isocurvature mode sources itself. By definition v\2(tk, tk) = 
and V22(tk,tk) = 1 at horizon crossing, since no interaction of the different modes has yet 
occurred (see also (3.2) and (3.4)). For the transfer functions of the adiabatic mode one finds 
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that vn = 1 and V21 = 0, since the curvature perturbation is conserved for purely adiabatic 
perturbations and adiabatic perturbations cannot source entropy perturbations. In order to 
better understand the role of the transfer functions, we can use the Fourier transformation 
of the perturbations [20] along with these last identities, to find 

/d 3 fc 
^ 3 /2 ^vima ] m (k)e lk ' x = Ci(i)Ofc) +«i2(Mfc)C2(l)(*Jfe)) 

/d 3 fc 
J^^7kV2 m al(k)e lk - x = v 22 (t,t k )( 2{1) (t k ), (3.12) 



where Cm(l) with m = 1,2 are the first-order adiabatic and isocurvature perturbation. 

Let us start by discussing the time evolution of v\2- Each one of the curves on the 
left-hand side of figure 3 corresponds to the time evolution of V12 for a different horizon- 
exit scale. At t = tk, i.e. when the relevant mode k exits the horizon, vi2k = since the 
isocurvature mode has not had time to affect the adiabatic one. Outside the horizon and well 
in the slow-roll regime of the sole dominance of the heavy field, isocurvature modes source 
the adiabatic one and the latter slowly increases. As time goes by, the heavy field rolls down 
its potential and the light field becomes more important. During this turning of the field 
trajectory, the mixing of the modes becomes maximal and violent isocurvature effects turn 
on. vi2k grows substantially during that era because of the increasing values of r] ± in (3.11) 
as well as the growing contribution of £>22fc> to become constant afterwards when the light 
field becomes dominant in an effectively single-field universe. 

Note that the earlier the mode exits the horizon the smaller is the final v\2k- This is 
opposite to the behaviour of the initial value, just after horizon-crossing, when the earlier the 
scale exits the horizon the more has its adiabatic mode been sourced by the isocurvature one 
at a given time t, and hence the larger is its i>i2fc- This can be understood by the evolution 
equations of v\2k and V22k in (3.11), showing that v\2k is sourced by V22k, which itself is a 
decreasing function of time, at least during eras when the universe is dominated by a single 
field (see the right-hand side of figure 3). If the equation (3.11) for vi2k did not depend 
on 77-*-, the vi2k curves would never cross each other since they would be similar and only 
boosted by their horizon-crossing time shift. It is the increasing value of rj^ or in other words 
the fact that the background is two-field dominated (and hence induces isocurvature effects) 
that results in the larger values of v±2k for larger k. 

On the right-hand side of figure 3 we show the evolution of 1)22- At horizon-crossing 
t = tk, V22k = 1> i-e. the isocurvature mode has not yet interacted with the adiabatic mode. 
Once outside the horizon and due to the dominance of the universe by a single field (but 
nonetheless an rj^ that is not quite zero), isocurvature modes are diluted while sourcing the 
adiabatic mode. During the turning of the fields isocurvature effects turn on, thus leading 
to an enhancement of ^22fc due to the varying value of x in (3.11). After the turning of the 
fields, the remnant isocurvature modes are again diluted to adiabatic modes and (for this 
model) at the end of inflation none are left. 

The transfer functions themselves may not provide much physical understanding of the 
system, since they are neither normalised nor can be recognised as any familiar physical quan- 
tities. Inspecting equation (3.12), one realizes that one can easily construct more physical 
quantities from the operators Cm(i) an d hence from v m2 , namely the ratios of the adiabatic 
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Figure 4. Left: The time evolution of the adiabatic R\ and the isocurvature i?2 ratios of power 
spectra as a function of the number of c- foldings N — iVgo f° r scales fcgo/10 (dashed line), Uqq (solid 
line) and fcgo x 10 (dotted line). Right: The time evolution of the spectral index (2.2) as a function of 
the number of e-foldings N — Nqq for scales k^o/lO (dashed line), k^o (solid line) and k$o x 10 (dotted 
line). Both plots are made for the model (3.9) with mass ratio m^/m a = 9. 



and isocurvature power spectrum to the total power spectrum: 

R (ClCl) _ 1 + (V12) 2 

1 (CiCi) + (C 2 C2> i + {v 12 f + {v 22 f 

D _ (C2C2) (^22) 2 ,„ 1 „s 

2 - (ClCl) + (C2C2) I + (y 12 f + (y 22 f { ' 

These are plotted on the left-hand side of figure 4 as a function of the number of e-foldings 
for different scales. One can clearly see that both ratios start as equal to 1/2 when the scale 
exits the horizon, while afterwards the adiabatic ratio R\ increases to reach 1 at the end of 
inflation and the isocurvature R 2 decreases to reach 0, for this particular model. During the 
turning of the fields violent isocurvature effects turn on and R 2 increases temporarily, while 
the adiabatic R\ has the opposite behaviour. 

On the right-hand side of figure 4 we plot the time evolution of the spectral index (2.2) 
of the power spectrum. The spectral index measures by construction the tilt of the power 
spectrum for different horizon-crossing scales and hence it depends on the horizon-crossing 
slow-roll parameters. For multiple-field models the power spectrum evolves during inflation 
even after horizon-crossing, and so does the spectral index. During the turning of the fields 
the spectral index increases, to remain constant afterwards. The earlier a scale exits the 
horizon the less negative is its spectral index ri£ — 1. This implies that the power spectrum 
itself decreases faster for larger horizon-crossing scales. This is due to the fact that except 
for the factor 1 + (v\ 2 k) 2 m the expression for the power spectrum there is also an inverse 
power of efc (see [20]). 

3.3 Spectral indices 

Finally let us discuss the scale dependence of the local /nl in terms of the relevant spectral 
indices. Equation (3.1) for an isosceles triangle implies that 

6 J_ f(k',k') + 2(^f(k',k) 

5 /nl { ^ 2?vdy)2 + 2{ v frmvdk y ^ ) 
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where 

f(k', k) = -2(k'kf 7 lnl v 12k 'Vi2k {g sr (k', k) + g iso (k', k) + g int (k', k)) . (3.15) 

For an arbitrary triangle configuration this is generalized as 

_6 _J_ fcf/(fci,fc 2 ) + perms. 

5 /NL (2n 2 ) 2 klV^k^ih) + perms.' [ ' 

The local /nl depends on a two- variable function /(fci,/^), with fei > &2- This is due to its 
super-horizon origin, which yields classical non-Gaussianity proportional to products of two 
power spectra. Hence one expects that the scale dependence of /nl can be expressed in terms 
of only two spectral indices, characterizing the function /. Notice that this is particular to 
the local case. In general the bispectrum cannot be split as a sum of two-variable functions 
and one anticipates that three spectral indices would be needed. 

The next issue to be resolved is which are the relevant spectral indices for /. The naive 
guess would be /(fei,^) = /(&i,o> fc2,o)(&i/&i,o) rifcl (&2/&2,o)™* 2 • We tested this parametriza- 
tion and we did not find good agreement with the exact value of /. Instead of that, we found 
that / is best approximated by keeping either the shape or the magnitude of the triangle 
constant. This statement can be expressed as 

/ K \ ™ K f \ ™ w 
f(ki,k 2 ) = f —) - , (3.17) 



where 

din/ , din/ 

n K = ts and n " = m — ( 3 - 18 ) 

dm A d in cj 

and 

U =k 2 = w^wr (3 - 19) 

The last equality is valid only for the isosceles case 7 = (see (2.10)). We dropped the —1 of 
the power spectrum spectral index definition to follow the definitions in [32, 33]. We added a 
tilde to indicate that these spectral indices are defined for the function /, not yet for the full 
/nl- In the next two sections we are going to examine the scale-dependence of /nl, changing 
the magnitude and the shape of the triangle separately, and verify assumption (3.17). 



4 Changing the magnitude of the triangle 

In this section we shall study the behaviour of /nl for triangles of the same shape but 
different size, see the left-hand side of figure 1. In figure 5 we plot the time evolution of 
/nl for equilateral triangles (the result would remain qualitatively the same for any isosceles 
triangle) of perimeter K = (3/2)&60 x 10 (top curve), K = (3/2)&60 (middle curve) and 
K = (3/2)/c6o/10 (bottom curve). The later the relevant scale exits the horizon the larger is 
its initial /nl as explained in the previous section, /nl grows during the turning of the fields 
due to isocurvature effects as described by (3.5) and (3.6), but by the end of inflation, when 
isocurvature modes vanish, it relaxes to a small, slow-roll suppressed value (see e.g. [19, 20]). 
In figure 6 we plot the final value of /nl (left) and the final value of the bispectrum (right) for 
equilateral triangles, varying K for values around K = (3/2)^60, within the Planck satellite's 
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Figure 5. The time evolution of /nl in terms of the number of e-foldings N — Nqq around the 
time of the turning of the fields, for equilateral (u> = 1) triangles with K = (3/2)fcgo (solid curve), 
K = (3/2)fc 6 o/10 (dashed curve) and K = (3/2)fc 60 x 10 (dotted curve), for the model (3.9) with mass 
ratio m<p/m a = 9. 
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Figure 6. Left: The relative change of the final value of /nl arbitrarily normalized to one at the 
smallest value of K on the figure, as a function of K for equilateral triangles (w = 1), calculated exactly 
(black curve), using the analytical approximation (4.1) (dashed red curve) and using the shape index 
(4.6) (dotted black curve). Right: The logarithm of the final value of the exact bispectrum, similarly 
normalized, as a function of K for equilateral triangles (uj = 1). Both figures are for the quadratic 
model (3.9) with mass ratio m^/m^ — 9. The vertical dashed line corresponds to K = (3/2)fc@o. 



resolution {k! jk ~ 1000). The later the scale exits the horizon, i.e. the larger K, the larger 
is the final value of /nl and of the bispectrum. 

The final value of /nl can be found analytically for the quadratic model within the 
slow-roll approximation. By the end of inflation V22,/ = so that (3.5) vanishes, while (3.6) 
can be further simplified to give some extra horizon-crossing terms and a new integral Qint 
that is identically zero for the quadratic potential (see (3.8) and [20]). For simplicity we give 
here the final value of /nl for equilateral triangles, 

/— \2 / ,|| . Vk \ , (- \3 / _L (^+vl-Xk)xk\ || 

3(ui2fe) [ek + Vk-Xk + ^j +(vi2k) Ufc- 1 v ± I + gfc + Vl 

hh,eq,f(k) = —2 . 

(l + (Vl2k) 2 ) 

(4.1) 

This formula is actually valid for any two-field model for which isocurvature modes vanish 
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Figure 7. Left: The time evolution of the conformal index uk (4.2) around the time of the turning 
of the fields, for triangles with w — 1 (black curves) and lj = 5/2 (red curves, below the u) = 1 curves), 
with perimeter K — (3/2)fcgo (solid curve), K = (3/2)fcgo/10 (dashed curve) and K = (3/2)fcgo X 10 
(dotted curve). The three points on the left correspond to the analytical values of the index as 
calculated from (4.7). Right: The final value of the conformal index uk (4.2) for triangles with oj = 1 
(black curve) and w = 5/2 (red dashed curve) as a function of K. Both figures arc for the model (3.9) 
with mass ratio m<t,/m a = 9. 



at the end of inflation and for which gi n t = 0, like for example equal-power sum models. 
Inspecting the various terms it turns out that although v\2k tends to decrease the value of 
/nl,/ as a function of k, it is the contribution of the horizon-crossing slow-roll parameters that 
wins and leads to an increase of the parameter of non-Gaussianity for larger horizon-crossing 
scales. Note that for equilateral triangles K is simply 3fc/2. 

We turn now to the spectral index nx- Using (3.18) with (2.11) and assuming that 
7 = and (3 = const., we can express uk in terms of the horizon-crossing time derivatives as 

s dlnf 1 din/ 1 /A 
n K (t;t kl ,t k2 ) = J + J (4.2) 

Cifci l — £fci ot k2 l — e k2 

Then /nl takes the form 

e i f(KA) (^y K{tk '° ,tk ' o) + 2^f(k' ,k ) U 

- t/nl - 



(4.3) 

Note that the ratios k' /k' Q = k/ko = K/Kq, since /3 = const. 

The above formula can be simplified in the limit of squeezed-triangle configurations, as 
well as in the equilateral limit. When one takes the squeezed limit a; 3 3> 1 (note that this 
would be true for (3 > 2/3), one finds: 

6 ! f(k' ,k ) / KY K H' tk ^~ ndt K?- ndt ^ )+2 

7/NL " 



6 / K yK(t k > o ,t ko ) 
/nl,o hH » ( 4 - 4 ) 



where 



(2^) 2 V c (k' )V c (k ) \K 

- 6 f ( K 

= -5 /NL '°Uo, 

n * . \ dm/ NL «91n/ NL 1 ain/ NL 1 

n K {t;t k >,t k ) = = — 1 — . (4.5) 

d m K at k i 1 — ey ot k 1 — e k 
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For the equilateral case to = 1, (4.3) becomes 



-r/NL 



i f{k' ,k' ) / ^yxH\ ) - 2nci \ 




(4.6) 



The conformal spectral index nx measures the change of /nl due to the overall size 
of the triangle, namely due to a conformal transformation of the triangle. For an isosceles 
triangle this is conceptually sketched on the left-hand side of figure 1, but it can be generalized 
for any shape, tir coincides with the ny NL of [32, 33] and grossly speaking it describes the tilt 
of /nl due to the pure evolution of the inflationary background (note that for an equilateral 
triangle this statement would be exact). 

On the left-hand side of figure 7 we plot the time evolution of the conformal spectral 
index for an equilateral u = 1 and an isosceles u) = 5/2 triangle that exited the horizon 
at three different times, namely for K = (3/2)^60 (solid curve), K = (3/2)&6o/10 (dashed 
curve) and K = (3/2) A^o x 10 (dotted curve). We plot the u) = 5/2 case only to demonstrate 
that the results remain qualitatively the same; we shall study the effect of different triangle 
shapes in the next section. The characteristic peaks that nx exhibits during the turning of 
the fields are inherited from the behaviour of /nl at that time and it is a new feature that is 
absent in the time evolution of the spectral index — 1 (see the right-hand side of figure 4) . 

In the context of the long- wavelength formalism we are restricted to work with the slow- 
roll approximation at horizon exit, so that the slow-roll parameters at that time should be 
small and vary just a little. This should be reflected in the initial value of the spectral index, 
which should be O(efc'). In another paper we will study models that do not necessarily satisfy 
this constraint, using the exact cubic action derived in [39] instead of the long- wavelength 
formalism. The earlier the scale exits, e.g. the dashed curve, the smaller are the slow-roll 
parameters evaluated at horizon crossing and hence the smaller is the initial nx- Indeed, 
using the definition (4.5) with (3.10) for the initial value of /NLm> we find for equilateral 
triangles 



which confirms the above statement. 

We notice that the initial, horizon-crossing, differences between the values of nx for 
the different horizon-crossing scales mostly disappear by the end of inflation, after peaking 
during the turning of the fields. The final value of the spectral index is plotted on the right- 
hand side of figure 7 and is smaller than its initial value. It exhibits a small running of 
0(10%) within the range of scales studied, inherited from the initial dispersion of its values 
at horizon-crossing. To verify that nx describes well the behaviour of /nl> we have plotted 
the approximation (4.6) in figure 6 where it can be compared with the exact result. We have 
also verified this for other inflationary models, including the potential 



studied in [20] , able to produce /nl of 0(1). The final value of the spectral index is two orders 
of magnitude smaller than the value for the quadratic model. This is related essentially to 
the fact that for the potential (4.8) the turning of the fields, and hence the slow-roll breaking, 



nK,in = 



(4.7) 



ek + vl 



W = b - b 2 (J 2 + 6 4 cj 4 + a 2 <j> 2 , 



(4.8) 
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Figure 8. The time evolution of /nl as a function of the number of e- foldings N — Nqq around the 
time of the turning of the fields, for a triangle with u> = 1 (solid curve), uj = 5/2 (dotted curve) and 
uj = 1000 (dashed curve), all with fixed perimeter K = 3/2kgo, for the model (3.9) with mass ratio 
m^/m a = 9. 

occurs near the end of inflation. This means that at the horizon-crossing times of the scales 
of the triangle, slow-roll parameters change very slowly and as a result the initial variation 
of /nl is much smaller than the one for the quadratic potential. As a consequence, the final 
tilt of /nl will be smaller. 

In Appendix A one can find the analytical expression for the final value of the spectral 
index hk, found by differentiating (4.1), but since the result is rather lengthy and does not 
give any further physical intuition, we chose not to present it here. 

5 Changing the shape of the triangle 

After studying triangles with the same shape but varying size in the previous section, we 
now turn to the scale dependence of /nl for triangles of the same perimeter but different 
shape, see the right-hand side of figure 1. In figure 8 we plot the time evolution of /nl during 
inflation for an equilateral uj = 1 (solid curve), an isosceles uj = 5/2 (dotted curve) and a 
squeezed uj = 1000 (dashed curve) triangle, all of perimeter K = (3/2) /ego, as a function 
of the number of e-foldings. The profile of the time evolution of /nl was discussed in the 
previous section. Here we are interested in the shape dependence of /nl- 

Although it is during the peak that the variation of /nl for different shapes is more 
prominent, its final value is also affected. On the left-hand side of figure 9 we plot the value 
of /nl at the end of inflation for triangles of perimeter K = (3/2)A;go, normalised by its 
value for the equilateral case (uj = 1), as a function of uj. The deviation of the values is 
small since it is related to horizon-exit slow-roll suppressed quantities. Within the long- 
wavelength formalism (or the 5N formalism) slow roll at horizon crossing is a requirement. 
Nevertheless, the important conclusion here is that /nl decreases when the triangle becomes 
more squeezed. This can be attributed to the fact that the more squeezed is the triangle, the 
more the fluctuation is frozen and behaves as part of the background when scale k' crosses 
the horizon. As a result the correlation between k and k' becomes less and the resulting 
non-Gaussianity is smaller (see also the discussion below equation (5.1).) 

An analytical formula can be found when applying the slow-roll approximation to ex- 
pression (3.1) at the end of inflation, when isocurvature modes have vanished. We perform 
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Figure 9. Left: The final value of /nl normalised by the equilateral /nl as a function of ui for 
triangles with K = (3/2)^60 calculated exactly (black curve), using the analytical approximation 
(5.1) (dashed red curve) and using the shape index (5.4) (dotted black curve). Right: The logarithm 
of the final value of the exact bispectrum normalised by the equilateral bispectrum as a function of 
uj for triangles with K = (3/2)^60 (black curve) and the same assuming /nl scale independent (red 
dashed curve). Both figures are for the quadratic model (3.9) with mass ratio m^/rria = 9. 



an integration by parts in the integral (see (3.8) and [20]; as before (ji n t = 0). More precisely, 
assuming that we are really in the squeezed limit k <^ k', the ratio J^/^k' becomes very 
large and we can ignore the equilateral terms that depend only on k' and not also on k. We 
also assume that the decaying mode has vanished to simplify the expressions for the Green's 
functions (see [20] and the discussion in section 3). Gyik'k can be set to zero as one can see in 
figure 3 (since it is basically equal to vyi and only involves times at the very left-hand side of 
the figure). Moreover, the same figure shows that V\2k/vi2k' ~ 1 (in the formula these ratios 
are always multiplied by slow-roll parameters, so that the deviation from 1 would be like a 
second-order effect), so that we find in the end 



fNL,sq,f — G22k'kfNL,eq,f(k') + 



1 - G 



22k'k 



1 + (Vl2k')' 



tk> + rf k ' + 



M. 

V\2k' 



Xk' (Vl2k' 



,(5.1) 



where /NL,e<?,/ is given in equation (4.1). The only quantity in the above expression that 
depends on the shape of the triangle is G^fc'fc, so it must be G22k'k that is responsible for 
the decreasing behaviour of /nl,sij,/- Indeed, increasing lo for a constant perimeter K of the 
triangle means increasing the interval i& — ty and hence decreasing the value of G^fc'fc (see 
the right-hand side of figure 3, since in the slow-roll regime V22 = £22)- This means that the 
interaction of the two modes becomes less important. In the complete absence of isocurvature 
modes G^fc'fc = and /NL,sg,/ takes its minimal value. It is only the isocurvature mode that 
interacts with itself and the greater is the difference between the two momenta the less is the 
interaction. Notice that the single- field limit of this result would correspond to G22k'k = 
and v\2 = 0. 5 

The decrease of /nl for more squeezed triangles seems contradictory to the well-known 
fact that the local bispectrum is maximized for squeezed configurations. In order to clarify 



inspecting equation (5.1) we notice that we do not recover the single-field squeezed limit result of [6]. 
This is to be expected since our /nl is the local one produced outside the horizon. As a consequence, we have 
only used the first and second-order horizon-crossing contributions coming from the redefinitions in the cubic 
action (see [20, 39]) as initial sources of /nl in the context of the long- wavelength formalism. The result of 
[6] on the other hand comes from the interaction terms in the Langrangian and hence is not the same. 
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Figure 10. Left: The time evolution of the shape index n u (5.5) around the time of the turning 
of the fields, for constant K = (3/2)keo (black curves) and K = (3/2)/c6o/10 (red curves, above the 
K = (3/2)kso curves), and for the shapes uj = 5/2 (solid curve) and ui = 1000 (dashed curve). The 
points on the left correspond to the analytical values of the index as calculated from (5.6). Right: 
The final value of the shape index n u (5.5) for constant K = (3/2)kgo (black) and K = (3/2)/c6o/10 
(dashed red) as a function of w. Both figures are for the model (3.9) with mass ratio to^/ttv = 9. 



this subtle point, we stress that the left-hand side of figure 9 is essentially the ratio of the 
exact bispectrum to the bispectrum assuming /nl as a constant (2.9) and hence the products 
of the power spectrum cancel out. We also plot on the right-hand side of figure 9 the final 
value of the bispectrum (2.9), normalised by the value of the equilateral bispectrum for 
triangles with K = 3&6o/2. Although /nl is maximal for the equilateral case, the bispectrum 
has the opposite behaviour, since it is dominated by the contribution of the products of the 
power spectrum, which leads to an increased bispectrum for the more squeezed shape. At 
the same time though we show that there is a small contribution of /nl itself, leading to 
smaller values of the bispectrum when compared to a bispectrum where /nl is assumed to 
be constant. 

In order to quantify the above results, we examine the shape index h u (3.18), assuming 
K = const and 7 = 0, 

~ = d\nf 1 1 8 In/ 1 2uj 

nuJ dty 1 - e k > 1 + 2u dt k 1 - e k 1 + 2co ' 1 ' ) 



In terms of n^, /nl takes the form 

6 1 f(k^k' ) + 2^f{kM{^ 

:/NL 



(5.3) 

where k'/k' = (1 + /3)/(l + (3 ) « 2w/(w + i) and k/k = (1 - /3)/(l - /3 ) oc l/(w + |). This 
can be further simplified in the squeezed region w > 1 to find 



7/NL 



f{k' M) 



nu,(t k , Q ,t ko ) 



(- 2 )V <(ti)n ,, 0)(SS i.)"<^>-v™ + ^»<-»»-' 



^ \ fkj(t k , o ,t ho )+n ( (t ko )-l g / ^ \ n u (t k , Q ,t kQ ) 

= ~hh,0 ( — = -r/NL,0 — , (5.4) 
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with 



^ _ dln/ NL _ dln/ NL 1 1 d In / NL 1 2uj 
d In uj dtf-i 1 — ey 1 + 2uj dtk 1 — 1 + 2a; 



The shape index n u describes the change of /nl due to the relative size of the two scales, 
namely due to how squeezed the triangle is, while keeping K constant (see the right-hand 
side of figure 1). 

We studied different squeezed triangle configurations with constant K, varying uj from 
u) = 5/2 to uj = 1000. On the left-hand side of figure 10 we plot the time evolution of the 
shape index. The negative values of the index signify the decrease of /nl as expected. As 
one can see from the figure, for the more squeezed triangle {uj = 1000) the initial value of n u 
seems to depend solely on the shape of the triangle and not on its magnitude, and even for 
the less squeezed triangle {uj = 5/2) the initial dependence on K is negligible. We can find 
the analytical initial value of n u by differentiating (3.10): 

1 4w {yj,) 2 

n>u),in = -, . n riK,in + -, . <J22fc'fc IT- {0.0) 

l + luj I + Iuj ek , + ^ 

For G22k'k = 0, which corresponds to the squeezed limit, n U) i n is proportional to the ini- 
tial shape index for equilateral triangles times a factor depending on the shape, which also 
becomes very small in the squeezed limit. 

Super-horizon effects, and especially the turning of the fields, result in a separation of 
the curves of n u of the same shape for different values of K, due to the dependence of the 
evolution of V\2k' on the scale k' . The turning of the fields increases the absolute value of n w , 
which is the opposite of the behaviour of the conformal index nx- The shape index depends 
on the transfer function Vi2k' (see Appendix A for an analytical approximation). The smaller 
K, the less does the final value of v\2k' change with respect to its initial value (see figure 
3) and hence the less the shape index is affected. Notice that the slow-roll parameters at 
horizon-crossing have the opposite behaviour: the smaller K, the smaller they are. Even 
though riu also depends on the slow-roll parameters, it is Vi2k' that most affects its evolution. 

On the right-hand side of figure 10 we plot the value of the shape index at the end of 
inflation. It exhibits a running of about 20% within the range of scales studied, somewhat 
larger than the conformal index. We have analytically computed the shape spectral index 
for models with cn n t = and with final gi so f = and give the result in Appendix A. 

The dotted curve in the plot on the left-hand side of figure 9 shows the final value of 
/nl approximated as a simple power law according to (5.4). Within the range of validity 
of our approximation uj^ 3> 1 it describes the exact result very well. We have also studied 
the shape spectral index for the potential (4.8). Similarly to n,Ki its value is two orders of 
magnitude smaller than the value for the quadratic potential, but the parametrization of /nl 
in terms of the shape index is in good agreement with the exact result for a larger range of 
uj > 3/2. 



6 Conclusions 

In this paper we studied the scale dependence of the local non-Gaussianity parameter /nl for 
two-field inflationary models. Multiple-field models with standard kinetic terms do not ex- 
hibit the strong scale dependence inherent in models that produce equilateral non-Gaussianity 
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at horizon-crossing through quantum mechanical effects. Nevertheless they are not scale in- 
dependent in general and the interesting question is whether we can profit from their scale 
dependence in order to observationally acquire more information about inflation. 

We have calculated /nl using the long-wavelength formalism. This constrains us to 
assume slow roll at horizon-crossing and hence the relevant quantities at that time should 
not vary much, including the scale dependence of /nl for any triangle shape. Indeed we 
confirmed that, by introducing the conformal spectral index uk that measures the tilt of 
/nl for triangles of the same shape but different size {K is a variable proportional to the 
perimeter of the momentum triangle). For the quadratic model with mass ratio m^/m CT = 9 
we find nx — 0.018, pointing to an almost scale-invariant /nl- 

We also studied the scale dependence of /nl while varying the shape of the triangle and 
keeping its perimeter constant, /nl exhibits the opposite behaviour of the full bispectrum, i.e. 
it decreases the more squeezed the triangle is (the momentum dependence of the bispectrum 
is dominated by that of the products of power spectra, not by that of /nl)- This variation 
is not related to horizon-crossing quantities, but rather to the fact that the more squeezed 
the isosceles triangle under study, the smaller the correlation of its two scales. We quantified 
this effect by introducing the shape spectral index n u , which for the quadratic model with 
r n<f>/' m (T — 9 is n u ~ —0.03 and has a running of about 20% (uj is defined as the ratio of the 
two sides of an isosceles momentum triangle). 

All our calculations have been done numerically in the exact background, assuming slow 
roll only at horizon-crossing, not afterwards. Nevertheless, semi-analytical expressions can 
be easily produced by directly differentiating /nl- If we do assume slow roll we showed that 
we can even simplify these expressions further and find analytical formulas for the final value 
of /nl and its spectral indices uk and n^, if the integral in /nl and the isocurvature modes 
vanish by the end of inflation, which is the case for example for any equal-power sum model. 

We used the two-field quadratic potential in our numerical calculations. This potential 
is easy to examine and allows for simplifications in the relevant expressions. Although its 
final non-Gaussianity is small, O(ejt), its general behaviour should not be different from 
other multiple- field inflationary models with standard kinetic terms, in the sense that the 
scale dependence of /nl should always depend on horizon-exit quantities and the evolution 
of the transfer functions during the turning of the fields. Indeed we have checked that for the 
potential (4.8) studied in [20], able to produce /nl ~ 0(1), the results remain qualitatively 
the same, although the values of the spectral indices are smaller due to the very slow evolution 
of the background at horizon-crossing in that model. 

Although the effect of the magnitude of the triangle on /nl had been considered before, 
analytical and numerical estimates were not available before this paper. In addition, it is the 
first time that the dependence of /nl itself (instead of the power spectra in the bispectrum) 
on the shape of the momentum triangle is studied. Using the long-wavelength formalism we 
have managed to study the two different sources of momentum dependence, i.e. the slow-roll 
parameters at horizon crossing and the evolution of the transfer functions, and to understand 
the role of each for the two different triangle deformations that we have studied. In summary, 
the later a momentum mode exits the horizon, the larger the slow-roll parameters are at that 
time and the larger /nl tends to be. In contrast, the final value of Vi2k' and the initial value 
of G^fc'fc, the two transfer functions that are the most important for /nl, are smaller the 
later the scale exits, which results in decreasing values of /nl- These two opposite effects 
manifest themselves in the two different deformations we have studied. When keeping the 
shape of the triangle constant and varying its size, it is the slow-roll parameters at horizon 
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crossing that play the major role in /nl and result in an increasing /nl for larger K. When 
changing the shape of /nl ; it is the correlation between the isocurvature mode at different 
scales, G 22k > k , that has the most important role, resulting in decreasing values of /nl when 
squeezing the triangle (i.e. increasing u). 

We have verified that the spectral indices of /nl {nx and n w ), which we introduced 
to describe the effect of the two types of deformations of the momentum triangle, provide 
a good approximation over a wide range of values of the relevant scales. In the models we 
studied their values are too small to be detected by Planck, given that /nl itself cannot be 
too big (or it would have been detected by WMAP). Models that break slow roll at horizon 
crossing could in principle have larger spectral indices, allowing one to detect them or rule 
them out using Planck data, but in order to study such models one would need to go beyond 
the long-wavelength formalism. We intend to study such models in a future paper using the 
exact cubic action derived in [39]. 



A Analytical expressions for the spectral indices 



By differentiating (4.1) and using (4.2) we can find the final value of nx for equilateral 
triangles in the slow-roll approximation, assuming that isocurvature modes have vanished for 
an equal-power sum potential (for which the gi n t contribution is zero, see (3.8) and [20]): 
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.Vl2k(vi2kXk - 2Vk) 



2e 2 k -3e k r 1 l + (vi) 2 + 5(rik L ) 2 -Q 



1 + (vuk) 2 /NL, eg ,/(l + (vuk) 2 ) 2 

+3v l2k l Vk L (3e k + 6r/| - 5 Xk ) - & J + 3(v 12k ) 2 ( W 22lk + 4(^) 2 - 2(e k + r,\- X k)(e k + 2 Xk ) 



<Xk{3e z k - 2W 221k + Ae kV l + 3(r,lf - 8fa£) a + V ll Xk - 3x1) + + vl ~ Xk) 



Ik 



+vl(el ~ (V l b 2 ) + W 221k (e k + 4) + ( V i)\2e k + 5^) - £ U + ^ + ^ Xk)Xk 



(A.l) 



where W 22 \ = (\/2e / k)W 22 \ / '(3H 2 ) . We have checked this approximation and we find good 
agreement with the exact conformal index for equilateral triangles. 

We repeat the calculation for the shape index n w (5.5), differentiating the squeezed /nl 
(5.1). Where needed we use the slow-roll approximation G^ 2k / k = —Xk'G 22k i k and G 2 ^ k i k = 
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G 22k ' k /2>. The result is: 

n uj,sqJ 

1 1 J 2G 22k i k V\ 2k i /, / . /„ . \ — x Xfe' 

1 m i F ^2^2 1 , o ^ i i F ^2 4 ~ ( w + ( 2 + w Vl2fc ' — l 



(A.2) 



+ «i2As' ( 3(e fe / + r/jl,) - 2xfc' ) + (v 12 k>) Vk 



-G22k'k 
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X k '(Sel-2W, 



22XW 



+4e fc ,7 ? [i, + 3(r7|!,) 2 - 7(^) 2 - 2j& - e vX w) + + vl - Xk>) + vl(4> ~ ^l?) 



+W 22W (e k , + rjl) + {rii,f{2e kl + 5^,) [Vk' + 



_L / _L , ( £ fe' + ^fc/ + XfeOXfe' 



+(vi2fc') 2 (2^22ifc' - 64 + (??f') 2 + 9(r/^) 2 " d' " 9 d'Xfc' + 8x1 ~ efc'^l + 5 X k> 



1 + (^12fc') 2 



K') + Hv k ') ~ 2e k > ~ 3(-k>%> ~ ~ 2{v 12k >) U k > + Vk'( e k' ~ 2%, - $Xk 



+{v\ 2 k') A \ W 22 ik' + (<*' - vl'Wl' + 3 (^) 2 + d' + 2 Xfc'( £ fc' + Xk')) ~ 2v 12k ,(£v + T}#(5e kl 
+2r)l + 3 Xk ')) + (vi 2k ') 2 (W 221k , + 2 (-e v {ey + r/J,) - 5(?^) 2 + (3e fc / + 377JJ, + X fc')Xfc- 



We have checked this approximation and we find good agreement with the exact shape index 
for ui > 3. 
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